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ABSTRACT 

We describe the effects of weak gravitational leasing by cosmological large scale structure on the 
diffuse emission of 21 centimeter radiation from neutral hydrogen at high redshifts during the era 
of reionization. The ability to observe radial information through the frequency, and thus three- 
dimensional regions of the background radiation at different redshifts, suggests that 21 cm studies may 
provide a useful context for studying weak leasing effects. We focus on the gravitational leasing effects 
on both the angular power spectra and the intrinsic, three-dimensional power spectra. We present a 
new approach for calculating the weak leasing signature based on integrating differential Fourier-space 
shells of the deflection field and approximating the magnification matrix. This reduces the problem of 
calculating the effect and higher order corrections to solving coupled systems of differential equations. 
This method is applied to reionization models of the 21 cm background spectra up to small angular 
scales over a range in redshift. The effect on the angular power spectrum is typically less than one 
percent on small angular scales, and very small on scales corresponding to the feature imprinted by 
reionization bubbles, due to the near-scale invariance of the angular power spectrum of the 21 cm signal 
on these scales. We describe the expected effect of weak lensing on three-dimensional 21 cm power 
spectra, and show that lensing creates aspherical perturbations to the intrinsic power spectrum which 
depend on the polar angle of the wavevector. The effect on the three-dimensional power spectrum 
is less than a percent on scales k < 0.1 h Mpc“^, but can be > 1% for highly inclined modes for 
fc > 1 h Mpc“^. The angular variation of the lensing effect on these scales is well described by a 
quartic polynomial in the cosine of the polar angle. The detection of the gravitational lensing effects 
on 21 cm power spectra will require very sensitive, high resolution observations by future low-frequency 
radio arrays. 

Subject headings: cosmology: theory — diffuse radiation — gravitational lensing 


1. INTRODUCTION 

After recombination at the surface of last scattering 
of the cosmic microwave background (CMB) at z = 
1100, most of the baryons in the universe are in the 
form of neutral hydrogen. Soon after the first galaxies 
and stars formed, they began to reionize the intergalac- 
tic medium. Understanding the astrophysical processes 
through which this happened is a major challenge in 
modern cosmology. Observations of the Gunn-Peterson 
trough in high-redshift quasars suggest that the univers e 
was reionized by redshift z ~ 6 fe.g. iFan et~ZI |2002|1 . 
Meanwhile, the large optical depth to Thomson scatter¬ 
ing of GMB photons observed by WMAP provides an 
integral constraint on reioni zation, indicating th at reion¬ 
ization began before z > 14 (iSoergel et al.l200'^ . At face 
value, the evidence suggests that reionization was a long 
and complicated process; it is a topic of intense current 
research. 

Since neutral hydrogen radiates in the rest frame at 
a wavelength of 21 cm due to the “spin-flip” hyperfine 
transition, 21 cm tomography has been suggested as a 
very promising probe of the history of reionization. This 
radiation should be observable at earth at low radio fre¬ 
quencies after being cosmologically redshifted. This 21 
cm radiation from the reionization epoch should carry 
immensely rich information about the reionization pro- 
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cess and the formation of structure. Furthermore, since 
by scanning through frequency one effectively is scanning 
through different redshifts of 21 cm emission, one has 
the remarkable ability to probe the evolution of the neu¬ 
tral hydrogen density with time. Early work studying 
the potential of 21 cm observatio ns to probe the high - 
redshift universe wa s do ne bv. e.g. iScott _ fc ReesI l|1990ll : 
iMadau et alJ lll997ll and iTozzi et al.1 (1200011 . 

Recent research has focused on the expected statistical 
signal of the 21 cm intensity from th e epoch of reioniza- 
tion (' EoR.l llZ aldarriaga et alJ |2004 iMorales fc Hewitil 
2f)f)4l IRha,radwaj fc Alii 1299^ . It was proposed by 
Zaldarriaga et all l|2004l that one may study the angu¬ 
lar power spectra of brightness temperature fluctuations 
at different frequencies (and cosmic epochs) in order to 
probe the details of the reionization history. An alterna¬ 
tive approach emphasiz es the full three-dimensi onal na¬ 
ture of the fluctuations ijMorales fc Hewitill2004ll . Angu¬ 
lar observations of the 21 cm sky stepped through differ¬ 
ent frequencies can be combined to form an image cube 
by mapping frequency to line-of-sight distance, and such 
three-dimensional image cubes can be used to probe the 
statistical properties of the full three-dimensional distri¬ 
bution of HI fluctuations at a particular epoch during 
reionization. The simple three-dimensional structure of 
the intrinsic fluctuations can be harnessed to remove 
foregrounds with distinct symmetry properties. Data 
analysis tech niques applying this app r oach have bee n 
developed bv IMorales fc Hewittl ll2004ll: iMoralesI l|2005ll : 
iBowma.n et al .1 ii2nn,^ and IMorales et al.1 l|200fi|l . Detec- 
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tion of the 21 cm signal from high redshift is among 
the major goals of planned future low frequency radio 
arrays, including the Mileura Widefield Array (MWA), 
the Primeval Structure Telescope (PAST), the Low Fre- 
qunecy Array (LOFAR), and the Square Kilometer Array 
(SKA). 

The expected statistical signal detected by future 21 
cm observations will be of a complicated nature, re¬ 
sulting from a combination of cosmological sources and 
astrophysical processes. The intrinsic 21 cm signal is 
determined by fluctuations in density, neutral fraction, 
and the spin temperature describing the relative pop¬ 
ulations of HI spin states, which is affected by the 
bac kground radiation field. iFurlanetto et alJ i2004all 
and lFiirlanetto et al.l ll2004bll have shown how inhomoge- 
nous reionization will imprint a shoulder feature on the 
power spectra; hence the power spectra can be an ef¬ 
fective discriminant betw e en po ssible reionization sce¬ 
narios. IBharadwai fc AlH II2004H have pointed out the 
importance of the effect of peculiar velocities, which in¬ 
troduce redshift space distortions and are caused by the 
infall towards (o r outfl ow from) density perturbations. 
iBarkana fc Loebl lj2005H have described how these pe¬ 
culiar velocity fluctuations lead to anisotropies in the 
intrinsic three-dimensional power spectrum P{k) with 
characteristic angular dependences on the orientation of 
k. These anisotropies may be exploited to separate out 
the effects of the density fluctuations and astrophysi¬ 
cal processes. The effects of radiative processes through 
the spin temperature have re centl y been c onsid ered by 
iPritchard fc Furlanettol ll2005ll and lHirat^ ll2005ll . 

These sources of fluctuations together determine the 
intrinsic power spectrum, and hence the autocorrelation 
function, at any given epoch (or redshift). As the 21 
cm photons emitted by HI stream through the universe, 
their paths are perturbed stochastically by density inho¬ 
mogeneities, and their observed directions in the sky are 
altered. Hence, gravitational lensing by large scale struc¬ 
ture presents another source of fluctuations and modi¬ 
fies the intrinsic power spectra and correlation functions. 
These weak lensing modifications are best thought of as 
secondary fluctuations; they are determined by density 
fluctuations along the photon path between emission and 
observation. Furthermore, they act upon existing fluctu¬ 
ations: without intrinsic brightness fluctuations charac¬ 
terized by the physics of the EoR, lensing would have no 
effect. 

The effect of weak gravitational lensing on diffuse 
backgrounds has been studied extensively in the con¬ 
text of the CMB. The effects on the temperature 
and polariza ti on power spectra hav e been discussed i n 
Seljaid (|199^ iZaldarria ga fc SeliakI lll9 98ll: i Hd ll200f]ll : 
CoOT^ 11200411 and iChalhnoTX LewisI ||200 !tIi . Further¬ 
more, the methods of using weak lensing to recon¬ 
struct i ntervening mass distr i bution s has been stu died 
in. e.g. I Zaldarriaga fc Se liakI lll999ll : iHul ll2001albll and 
iHirata fc SeliakI l|2003ll . Future observations of 21 cm 
backgrounds will present a new source of diffuse emis¬ 
sion that will be subject to weak lensing by large scale 
structure. In contrast to the CMB, which presents a sin¬ 
gle, narrow surface at one redshift, the 21 cm emission 
from different redshifts provides multiple source planes 
that may be observed in order to study lensing effects. 
Indeed, the three-dimensional nature of 21 cm tomogra¬ 


phy provides a new and interesting context for lensing 
sudies. Som e initia l steps have been taken in this di¬ 
rection. iPenI ll2004ll has suggested a method to extract 
information on the projected lensing convergence from 
future 21 cm surveys. 

In this paper we focus on the expected weak grav¬ 
itational lensing effects on the angular and three- 
dimensional power spectra of 21 cm emission from the 
era of reionization. In m we motivate and derive a new 
differential method for calculating the lensing modifica¬ 
tions. Additionally the lensing effects in certain limiting 
cases are explored. In isi we apply these methods to 
models of the 21 cm angular power spectra. In 0 we 
discuss the lensing distortions to the three-dimensional 
power spectra of 21 cm fluctuations, and show that they 
produce asymmetric perturbations. These effects have 
not been considered in previous work. In (l^we compare 
the magnitude of these effects to expected sensitivites of 
a future SKA experiment. Finally, we conclude in 0 


2. DIFFERENTIAL EQUATION APPROACH TO WEAK 
LENSING 


The case of weak gravitational lensing effects on the 
CM B power spe c trum has been thoroughly explored 
fe.g. iSeliakI 119961: IZaldarriaga fc SeliakI 119981: IHuI 120001: 
iChallinor fc Lewial200,^ lCooravll2004jl . The 21 centime¬ 


ter background from the EoR presents a distinct case in 
which to study the weak lensing effect. The two cases dif¬ 
fer significantly in the distribution of fluctuation power 
on small scales versus large scales. In the case of the 
CMB, the intrinsic power is dominated by the acoustic 
peaks on scales of order I ~ 100, while small scale power 
I > 1000 is suppressed due to Silk damping. The con¬ 
tributions to the temperatures fluctuations are confined 
to angular scales larger than arcseconds. In the case of 
the 21 cm background fluctuations from EoR, the power 
spectrum is dominated by small scale power on scales 
I > 1000. In order to fully understand the effects of 
weak lensing by large scale structure, one must consider 
many more decades of scale, and much smaller scales, 
than in the case of the CMB. 

This simple observation has implications for how one 
calculates lensing effects on the 21 cm background. Ex¬ 
isting methods for calculating the lens effects on power 
spectra of diffuse backgrounds are similar in that they are 
perturbative in the deflection angle field, which maps the 
source plane into the observed image plane. For example, 
in the method of Hu (2000), and as extended to higher 
order in Cooray (2004), the lensed temperature field is 
expanded in the deflection angle as a Taylor series. 


f{e) = T{e + se{e)) 

« T{e) + V,T(0)^6»* + iVaV6T(6»)56i“(56i'' -L ... 

( 1 ) 


The deflection angles are coupled to higher and higher 
derivatives of the temperature field. In cases where there 
is dominant small scale power, and thus significant fluc¬ 
tuations with wavelengths smaller than the typical de¬ 
flection angle on arcminute scales, one expects that ex¬ 
pansions of the temperature field in which only a few 
terms are retained are insufficient to describe the lensing 
effect. A low-order perturbative expansion of the un- 
lensed temperature in the deflection angle is an inaccu- 
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rate approximation of the lensed temperature field if the 
temperature is not sufficiently smooth on scales smaller 
than and of order the typical deflection angle. In order 
to characterize the lensing effect, it is apparent that one 
must carefully take into account higher order terms in 
these methods. 

There is a caveat in this reasoning. Since ultimately 
we wish to understand the lensing effect on the power 
spectrum, the lensing effect comes through only the two- 
point function of the product of the temperature at one 
observed point and another point. One can show that 
the calculation depends only on the difference between 
the deflection angles at the two points, rather than the 
individual magnitudes of each deflection angle. This is 
simply a consequence of the fact that an overall, coherent 
deflection of the temperature pattern on the sky will not 
change the pattern. Hence the relevant physical quantity 
is the variation of the deflection field over the scale of 
temperature fluctuations. 

A series solution for the lensed power spectrum is 
obtained by forming the lensing correlation function 
{T{9)T{9')), and performing the transform analytically. 
The result is 


cr 



(fki P cPkn 
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X [fci . (i - fci - . . . - . . . [fc„ . (i _ fci - . . . - kr,)f 

‘■-I’* 

( 2 ) 


where R = (l/47r) f dlnkk^C^ llCooravll2004ll . In order 

to compute the result to the fVth order in Cf, one must 
expand the exponentials and retain all integrals with up 
to N factors of the lensing potential powe r sp e ctruin , C^. 

In th e harm onic space expansions of iHd ll200nll and 
iCooravl ll2004ll . this dependence on only the relative de¬ 
flection angle manifests itself in cancellations between 
positive and negative integrals at a given order in the 
lensing potential power spectrum. When there is signif¬ 
icant small scale power, the individual terms at a given 
order can become very large, but the total contribution 
of these terms may be very small due to these cancel¬ 
lations between positive and negative terms. A correct 
computation of higher order contributions to the lensed 
power spectrum must carefully account for delicate can¬ 
cellations between such terms. Hence, the presence of 
significant small-scale power alone does not guarantee 
that appropriately summed higher-order corrections will 
be important. The degree to which such cancellations oc¬ 
cur generally depends on the shape of the intrinsic power 
spectrum. Nevertheless, one must check these high or¬ 
der contributions in order to determine whether or not 
the lensing calculation has converged properly on small 
angular scales. 

In the harmonic space expansion, however, the calcula¬ 
tion of higher order terms becomes numerically difficult. 
Whereas the leading order lensing contribution (at order 
cf) is a simple integral over the two-dimensional space 
of the lensing potential wavevector, the second order cor¬ 
rection requires the evaluation of a four-dimensional inte¬ 
gral. In general, the computation of the {Cf‘^)^ correc¬ 
tions requires the evaluation of a 2A^-dimensional numer¬ 


ical integral over the unlensed temperature spectrum and 
N lensing potential power spectra. In order to compute 
the lensing effects on the 21 cm spectrum, in which there 
is significant small scale power in several decades, these 
integrations must be done over many decades of scale. 
At each order, one must correctly handle the delicate 
cancellations between these multi-dimensional integrals 
over decades in scale. This procedure for accounting for 
higher-order effects becomes numerically exhaustive and 
slow, and may lead to inaccurate estimates. 

The other approach used in calc ulating the le nsed 
CMB power spectrum is de scribed in iSeliald l|199fill and 
iZaldarriaga fc SeliakI lll998ll . In this formulation, the 
two-point autocorrelation function of the lensed temper¬ 
ature field is computed numerically, and the transform 
of this result is computed to calculate the power spec¬ 
trum. This is also essentially an expansion in the full de¬ 
flection angle coupled to derivatives of the temperature 
field. Indeed, one can show that these two approaches 
are equivalent when one does not terminate either se¬ 
ries or expands each to the same order in the deflection 
dispersion. A downside to this method is the need to 
transform back and forth between position and Fourier 
space. This involves numerical integrations with Bessel 
functions which become highly oscillatory at small scales. 
Furthermore, the transforms obscure the effects of the 
lensing power spectrum on modifying the temperature 
power spectrum. A purely harmonic calculation has the 
appeal of transparency, as all relevant quantities interact 
in the Fourier domain. 

The difficulties apparent in these existing methods mo¬ 
tivate us to develop a new method to handle the lens¬ 
ing calculation in the case of the 21 cm background. In 
this section we present a new method using a differential 
equation approach. This method avoids the pitfalls of 
the above methods. First, the differential approach al¬ 
lows a straightforward computation of higher order cor¬ 
rections without requiring evaluation of irreducible high¬ 
dimensional integrals. The differential approach allows 
us at once to sum up many terms from all orders in 
the series. Cancellations between terms, corresponding 
to differences in deflections, are carefully done shell by 
shell in Fourier space as a function of a single wavenum¬ 
ber cutoff A. Extensions of the calculation to include 
finer approximations are straightforward to derive. The 
systems of differential equations are solved entirely in 
harmonic space, avoiding highly oscillatory Bessel trans¬ 
forms back and forth between position and harmonic 
spaces. Furthermore, the harmonic space solution allows 
for illuminating analytic simplifications in certain limits. 
Our new approach is easy to implement, requiring only 
one-dimensional integration and a differential equation 
solver. 


2.1. Formalism 

Our analysis will be carried out in the formalism of the 
flat-sky approximation, which is sufficient on small an¬ 
gular scales and small sections of the sky. Since we will 
apply this method to the 21 cm emission, where most of 
the power is at small scales, this approximation should 
be sufficient for our purposes. We expect that correc¬ 
tions due to the curvature of the sky will be small in the 
regime studied in this work. In general, for high-precision 
and low multipoles, one should take into account curva- 
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ture effects (e.g. iHd (1200(1) . IChallinor fc LewisI (120051) '). 
In the ffat sky limit, an observeable X on the sky can 
be expanded in Fourier modes rather than spherical har¬ 
monics, e.g. 

’ ( 3 ) 

with power spectra 

{X{l)X{l')) = {27rf (4) 


The leasing modes are usually approximated as Gaus¬ 
sian random variables. The fluctuations of lens quanti¬ 
ties, e.g. the projected potential, deflection angle and 
effective convergence, are described by the power spec¬ 
tra, which are integrals over the matter fluctuations along 
the line of sight. The convergence power spectrum to a 
source at comoving radial distance and redshift Zg in 
the small scale limit using the Limb er approximation is 
(e.g. iBartelmann fc Schneidej [20^: 




dr 


W{r, TgY 
a?{r) 


Ps 



( 5 ) 


where W{r,rs) = x(t s — r)/x(.fs)■, Ps denotes the matter 
power spectrum, and y is the comoving angular diameter 
distance. The power spectrum of the projected lensing 
potential is related simply through the Laplace equation, 

2k( 6/) = V2<(>(6>), so that = ACr. _ 

use the code CMBFAST (iSeliak fc Zaldarriagal 
Il99(1) to generate the transfer functions and to compute 
the lensing convergence power spectrum. On small scales 
it is important to account for the non-linear evolution 
of the matter power spe ctrum. We use the method of 
iPeacock fc Doddil (Il99(1) to calculate the non-linear ap¬ 
proximation. The linear and non-linear power spectra 
for several redshifts in the range of interest are plotted 
in Figure n These will generally depend on the cosmo¬ 
logical model. For the rest of this paper, we adopt a 
ACDM model with these values of the cosmological pa¬ 
rameters: = 0.3, Ob = 0.046, Oa = 0.7, h = 0.70, 

with mass fluctuations on scales of 8 h Mpc~^, as = 0.88 . 
in accordance with the latest data (iSoergel et al.ll200'^ . 


2.2. Derivation 

In this section, we derive a new method for calculating 
the lensing effect on the power spectrnm. Our strategy 
differs fundamentally from previous approaches. Instead 
of finding a solution perturbative in the full deflection 
angle, we expand in a differential shell of modes of the 
deflection field. This leads us to a differential equation 
solution for the lensing effect in Fourier space. We show 
how this strategy leads to a hierarchy of systems of differ¬ 
ential equations. A method for truncating the hierachy 
in the limit of weak lensing is proposed, leading to simple 
closed systems as first and second approximations. We 
explore various limits of these equations to gain intuition 
of the lensing process in Fourier space. 

In what follows, we will adopt the following conven¬ 
tions. The unobserved, source-plane quantities are de¬ 
scribed by non-tilde notation, e.g. the unlensed, intrinsic 
temperature power spectrum is Cf and the source tem¬ 
perature is T{f3). Quantities that are affected through 
lensing are described by tilde notation. Hence, the lensed 



Fig. 1.— Angular power spectrum of the lensing convergence. 
The solid curves are the non-linear spectra for 2 : = 1100, 15, 12, 
and 8. The dashed curves are the corresponding linear spectra. 


temperature power spectrum is Cf and the lensed, ob¬ 
served temperature is T(6). 

The lensed temperature field is related to the unlensed 
field through the deflection angle. 

f{9) = T{i3) = T{e + 6e{e)). (6) 

The deflection field can be decomposed into Fourier com¬ 
ponents. 

mO) = I « ( 7 ) 

Now we introduce the notion of a harmonic-space “mo¬ 
mentum” cutoff on lensing quantities. Define the cutoff- 
dependent deflection field including lensing modes up to 
wavenumber |fc| < A, 

DeiA)= t ( 8 ) 

Jo ( 271 -)^ 

and the cutoff-dependent “source position” (3t^[0) = 6 + 
Dg{A). The temperature field lensed by lensing modes 
with wavenumber up to |fc| < A -|- SA is 

Ta+5a( 0) =T(0 + £)e(A-f (5A)) 

= T(/3a + de) ^ ’ 

where we have defined the differential deflection element 
including only lensing modes with A < |fc| < A -|- 6A: 

de = [De{A + SA)-DeiA)] (10) 

One may interpret the cutoff A as a low-pass filter on 
lensing quantities that allows one to calculate the lens¬ 
ing effects of only long-wavelength modes on the lensed 
quantities. In any numerical calculation, there is a prac¬ 
tical cutoff on the maximum frequency of a Fourier mode 
included in the calculation of lensed quantities. We have 
simply introduced this cutoff explicitly, and use it as a 
tool to regulate the calculation of the lensed quantities. 
We will show how the lensing problem reduces to a sys¬ 
tem of differential equations in the cutoff, and this system 
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will be evolved to a high enough A so that the calculation 
in the regime of interest has converged and thus is inde¬ 
pendent of higher frequency spatial modes of the lens. 

The lensed temperature correlation function, lensed by 
deflection modes fc < A is formed as: 


^A+5A (|A6/|) = {fA+SAi9)fA+SA{e')) 

= {T{(3a + d)T{(3'j, + d)) = {T{(3a)T{(3' + Ad)) 

1 . 


= (TAiO) 


T{I3'a)+ T„(/3;)Ad” + -r„(/3^,)Ad“Ad* 


= M) - i( J^T(/3A)^T(/3;))(Ad“Ad* 


( 11 ) 


where primed and unprimed quantities correspond to the 
coordinates 6 and 6 ', respectively, and A6 = 9—6' is the 
displacement vector. Gradients of the temperature field 
with respect to the source position (3a are denoted by 
subscripts, e.g. dT/d(3\ = Ti{(3A)- We have expanded 
to second-order in the relative differential deflection an¬ 
gle Ad = d' — d and used the translational invariance 
property of the average. The Einstein summation con¬ 
vention for repeated indices is implied in each term, and 
i,j = 1, 2 = 9x,dy. Here we have expressed the temper¬ 
ature lensed by lensing modes |fc| < A -|- ^A in relation 
to the lensed temperature just including lensing modes 
|fc| < A, both at the observed position, 9. In the last two 
terms, the first ensemble averages are over the temper¬ 
ature and lensing modes up to |fc| < A, and the second 
averages are over the lensing modes in the Fourier shell 
A < |fc| < A -I- <5A, and we have assumed the lensing 
mode amplitudes are distributed as Gaussian random 
variables. This arrangement of terms emphasizes that 
the lensing effect depends only on differences between 
deflection angles at different points, not on their abso¬ 
lute values. This is simply a consequence of the fact that 
the correlation functions only depend upon differences in 
position. Since we will divide through by SA and take its 
limit to zero, the neglected terms in the Taylor expansion 
involving higher order products of the relative deflection 
will vanish; hence, this expansion becomes exact. 

In order to evaluate the derivatives on the source plane 
with respect to /3 a we must make a change of coordinates 
to relate them to the observed coordinates 9: 


|fe| = A. The inverse is then simply, 

M=_ i _ A + «-7i -72 \ 

(1 + k )2 - (72 + .^ 2 ) -72 l + ^ + 

+ 72 \ 

72 K-^lJ ■ 

(15) 

With these substitutions, the equation for the correla¬ 
tion function expressed in observed coordinates is simply 


CI+,a(A0) - CliAd) 

1 


2T 


M(dj 


Ta 


Mid', 


f^){Ad'Ad‘' 


(16) 


This equation relates the lensed correlation function with 
lensing cutoff A -|- (5A to the lensed correlation function 
with cutoff A; i.e. it is the change from integrating over a 
differential shell of lensing modes. The correlations of the 
differential deflections will be proportional to i5A in the 
limit 6 A 0. We will see that this leads to a differential 
equation for the lensed correlation function as a function 
of the cutoff A, or for the lensed power spectrum when 
we transform directly to harmonic space. 

There is an alternative way to derive the above re¬ 
sults that motivates another intuitive picture. The full 
deflection angle tells us how to map the source plane 
temperature to the image plane temperature. One 
can imagine intermediate steps in this mapping, which 
are image planes that are lensed through the cutoff- 
truncated deflection field, I3e(A), each labelled by the 
cutoff wavenumber, A. Just as the lens equation relates 
the intensities between the source plane and the lensed 
image plane, one can derive a relation between these 
intermediate lensed image planes. Gonsider the image 
lensed by modes up to A-I-JA, Ta+sa{9). One can relate 
this to a previously lensed image plane, Ta{9), lensed 
only by lensing modes up to A, through the lens equa¬ 
tion. Let the vector Xg be the effective infinitesimal 
deflection angle that maps the image at A to the image 
at A -I- 5A. 


fA+sA{9) = fA{9AXe) (17) 


The effective deflection Xg can be determined by map¬ 
ping each side to the source plane: 



de^ d 
8(31 dOi 




( 12 ) 


The matrix M^^A) = 86 ^ /8(3]^ is the inverse of the 
magnification matrix of the A-filtered lensing quantities 
Al{A) = 881/80’': 


M/Al 


80^ 88 \ 
88 \ 90’’ 



The magnification matrix is explicitly: 


(13) 


A^A + ^.ii \ ^A + k + 7i 72 \ 

\ 4>,12 1 + 4',22 ) \ 72 l-|-K-7iy 

(14) 

where we have defined the usual convergence and shears, 
K = 5 (</',ii + 8 , 22 ), 71 = - 4 >, 22 ), and 72 = 8 , 12 , 

and all quantities are implicitly cut off at lensing modes 


fA+5A{9)=T{9 + Dg{A)+dg) (18) 

fA{6 + Xg)=T{9 + Xg + Dg+x{A)) (19) 

Since the two expressions are equal, Xg + Dg+x (A) = 
Dg{A) + dg. Next the deflection field is expanded in 
the infinitesimal vector X: Dg^^{A) — Dg = {8bDg)X’’. 
Then we have 

d)^ = {3/+ 8 bD^{A))X^ = A/{A)X^ (20) 

so that Xg = M(/d’^. Thus the effective incremental 
deflection between two images labelled by their lensing 
cutoffs, A-|-<5A and A, is simply the incremental deflection 
dg rotated and magnified by the inverse magnification 
matrix of the lensing modes up to A. This relation can 
then be used to derive the above Eq. m, by expanding 
the temperatures at each position, Ta+ 5 a(^)) in A! to 
relate it to the previous lensed image, Ta{9). 
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2.2.1. The M = I approximation 

As a first approximation, if the lens is sufficiently weak, 
we may approximate the inverse magnification matrix as 
identity M^(A) ~ 5^. In doing so, terms proportional to 
K, 7 i, and 72 from Eq. m are neglected. This is justi¬ 
fied in the weak lensing regime in which the elements of 
the distortion tensor are much smaller than unity. Note 
that this approximation depends only on the strength 
of the lensing field and not the underlying temperature 
power spectrum that is being lensed. We can see from 
Figure^that the variance in the convergence is much less 
than unity. In this approximation, the lensing equation 
then simplifies to: 

=-i(ftTA5;r;)(Ad*Ad“') (21) 

where di = d/d9^ denotes derivatives with respect to the 
observed coordinate. This is evaluated by decomposing 
each quantity into Fourier components and averaging the 
coefficients as (TA(i)*TA(Z')) = {2'k)'^C'[[K) 5'^{l — V) and 
{4>{k)*(j){k')) = (27r)^C'^i5^(fc — k'). The lensed power 
spectrum is obtained by applying a Fourier transform 
(/d2A0e-'' (»-«') ) 


ACr(A) = ( J d4>K [A • (Z - A)]^ (7|T_^| (A) 

-[A-Zf cr(A)) 

( 22 ) 

and taking the limit JA ^ 0. In the limit, this simpli¬ 
fies to an ordinary differential equation for evolving the 
temperature power spectrum with the cutoff A. 


dCfiK) 
(Zln A 


PA^Ci 


27T 


Cl (A) cos 


27T . 


'K 


— cos 


c;[_ai(a) 


(23) 


where |Z — Ap = -|- A^ — 2ZA cos((/)a)- This solution 

contains contributions from all orders in Cf, and the er¬ 
ror contribution begins at order 0{{C^)‘^) in the sense of 
the series solution, Eq. ©• This approximate differen¬ 
tial equation is simple to solve, and should be adequate 
in the limit of very weak lensing, when the magnification 
matrix does not deviate significantly from identity. 

There are two limiting cases of interest. For A <C Z, 
and sufficiently smooth power spectra, we may expand 
the right-hand side in x = A/I and perform the integral. 
In this limit, the two terms in the integral cancel, and 
the residual terms are suppressed by a factor of x'^: 


dCf{A) 
din A 


PA^Cj r 

47r 


Ox^ 




A^Cj 

An 




d\nP J 
i d?Cf{A) 

8 dlnZ2 


+ 0{x^) 
(24) 


We see that in this limit, in which the differential shell 
of lensing modes has a wavelength much longer than 


that of the temperature fluctuation mode, the leading 
order lensing effect, proportional to in the first line, 
vanishes. Physically, this corresponds to a picture in 
which the deflection field is coherent over the scale of 
the temperature fluctuations. Temperature fluctuations 
on small scales are gravitationally deflected together, re¬ 
sulting in an overall shift that does not change the fluc¬ 
tuation power on those scales. In this limiting case, 
the vanishing of the leading term is independent of the 
power in the lensing mode. Thus, even though the RMS 
amplitude of the incremental deflection angle may be 
large compared to the fluctuation with wavenumber Z, in 
this particular limit the leading lensing effect vanishes. 
The remaining fractional lensing effect is proportional to 
A®C'^/(87r) = A^C//{2n), which is the power per loga¬ 
rithmic interval in the convergence. As one can see from 
Figure ^ this is a small dimensionless fraction, which 
declines steeply on large scales as A ^ 0. This intu¬ 
ition demonstrates that low A lensing modes are rela¬ 
tively unimportant, unless I ~ A. The lensing behavior 
then will generally depend on the power of lensing fluc¬ 
tuations at that scale and the shape of the temperature 
spectrum. 

It has been pointed out bv iHul (1200(11 that, since the 

power in the deflection field A^C^ peaks at a low multi¬ 
pole, A ^ 40 in ACDM, the lensing effect on even high-/ 
temperature modes must arise from large-scale power in 
the deflections. This reasoning, however, fails to fully 
account for the cancellations among deflection angles 
over the scale of the temperature fluctuations. We have 
shown, in Eq. lEH), that in this limit, A ^ Z, the rele¬ 
vant quantity determining the differential lensing effect, 
after accounting for coherent shifts, is actually the power 
per log interval of the convergence, A®C'^/(87r), which is 
quite small on large scales and peaks at much smaller 
scales (Fig. We will see in (j^that, in the case of 
21 cm fluctuations, the contribution by large-scale lens 
modes (A < 50) to the total lensing effect on smaller 
scale temperature fluctuations is negligible. 

In the opposite limit, A ^ Z, one may expand the right- 
hand side in y = Z/A. The result in terms of the variance 
per logarithmic interval, A^^ = /{2,n), is 


c/Ar(A) A^ci 


cZln A 
Ay 


2n 


-y 


Af(A)--Ar(A) 


,2/ 1a2t^a'i I ^d'^Al^{A) 
--Aa (A) +4 


fc=A 


0{y^) 


(25) 


In this limit, there is a small wavelength deflection mode 
lensing a long wavelength temperature mode. Power on 
the scale Z is modified, and the magnitude of the ef¬ 
fect is determined by a combination of the convergence 
power at the scale A and the temperature spectrum. The 
overall effect is determined by the convergence power, 
A^C^ ~ A^Ca, suppressed by a factor y^; hence we ex¬ 
pect that the lensing effect will shut off in the limit, 
so long as the convergence power is slowly varying on 
sufficiently small scales. From Figure ^ we see that, 
whereas the convergence power is increasing on large 
scales A <C 10^, it is fairly flat or decreasing on small 
scales A > 10^. The leading order effect is determined 
by the distribution of temperature fluctuations on large 
scales (Z) versus small scales (A). If the fluctuations are 
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dominated by small scale power, the small 

scale lensing modes produce power on scale 1. In this pic¬ 
ture, a short wavelength lensing mode can interact with 
a temperature mode of the same wavelength to produce 
long wavelength power. If there is negligible small scale 
power, <C A;^^, small scale lensing modes remove 
power in fluctuations at scale 1. This is intuitively the 
case, for example, for the damping tail of the CMB, in 
which lensing deflections create smal l scale fluctuation s 
against a long wavelength gradient llZaldarriagal l200n(l . 
In this case, there is little intrinsic power in small scale 
temperature modes I > 4000, and small scale (A) deflec¬ 
tions of large-scale gradients produce small scale temper¬ 
ature fluctuations on scales I ~ A. Consequently, since 
lensing preserves total variance, power is removed from 
the acoustic regime I < 2000. The above limit demon¬ 
strates this overall shifting of fluctuation power. 

In these limits we were able to replace an integral over 
Cf (A) with its local derivatives. Hence these formulae 
will be useful in a numerical implementation for suffi¬ 
ciently smooth power spectra. Furthermore, these local 
derivatives of the temperature power spectrum with re¬ 
spect to I emphasize the importance of the shape of the 
power spectrum in determining the lensing effect. 

We note that in cases in which Eq. PI is a sufh- 
cent approximation, and the lensing modifications to the 
power spectrum are small for all A, one may make the 
replacement Cf (A) Cf on the right-hand side of (1^ . 
Then the right-hand side can be immediately integrated, 
and the solution is 

Cf[K)-Cf= t^ct{-{k-lfcr 

which, as A ^ oo . is the first-order harmonic-space solu¬ 
tion of 113 (E^, or, equivalently, the first order terms 
of the full series solution, Eq. obtained by the sim¬ 
plest expansion of the temperature field to second order 
in (full) deflection angles. Hence, we see that the solu¬ 
tion to Eq. 113) includes this most basic approximation, 
but sums up higher order contributions by having the 
current lensed power spectrum Cf (A) on the right hand 
side of Eq. (|13. 

The difference between the lensing solutions, Eq. 113 
and Eq. PI is in the nature of the approximations. 
In Eq. it^ the terms that are dropped to make the 
approximation are those corresponding to the higher or¬ 
der terms in the Taylor expansion, Eq. o, of the form 

.. .dif^T{6)59'^^ .. .59^^. In Eq. (1^ . the terms 
that are dropped are proportional to the elements of the 
distortion tensor, and are of order 0 (k(A)) smaller com¬ 
pared to the terms retained, i.e. those corresponding to 
the identity part of the inverse magnification matrix. It 
is in this sense that Eq. PI sums the series Eq. © in 
a different way than Eq. Pi. 

2.2.2. Poisson Power Spectrum 

The solution to Eq. PI is simple in the special case 
of a Poisson power spectrum. In the Poisson case, Cf = 
constant is independent of 1. Lensing acts to remap un¬ 
correlated points, hence the lensed power spectrum re¬ 
mains Poisson. One can then take Cf (A) as indepen¬ 


dent of I and remove it from the integral in Eq. P|. 
The solution to the equation is then: 

C[ = Cf exp( J ^4Cfe") « Cfil + Ainf) (27) 

where is the power spectrum of the convergence. 
Thus the temperature spectrum is modified by multi¬ 
plication by a constant factor determined by the total 
power in the convergence. However, the spectrum re¬ 
mains Poisson, i.e. independent of I, under the action of 
lensing. 

We can understand this by considering a simple model. 
Consider a region with temperature fluctuations, de¬ 
scribed by a Poisson spectrum up to some minimum scale 
^ magnified by a constant convergence so that 

T{9) = T{p,~^9), where p is the magnification factor. 
Then the lensed correlation function is 


{f{9^)f{9^)) 


{T{p-^9^)Tipi-^9f) 

d^k 

J 

pf^raa.x/ fd' (]j‘^ I 


(27r)- 


2^T ilAe 


(28) 


In the case of a Poisson spectrum, then Cf = const x 
cf on scales I < kmax/t^- One can see that the small 
scale cutoff kmax is rescaled in order to preserve the total 
variance. 


2.2.3. Scale-Invariant Power Spectrum 

Suppose at some cutoff A, the variance Af^{A) = 
PCf{K)/{2'K) is approximately constant independent of 
wavelength, at least within the range of integration. 
Equation E3 simplifies as: 

dAr(A) A6cl,2T..^^0 ifA<Z, 

din A - 4n | ^ [l _ (^)2] if A > Z. 

(29) 

This shows that small scale temperature modes are not 
affected by lensing by larger scale lensing modes at scales 
A < /, independently of I in this regime. Furthermore, 
much shorter wavelength lensing modes A ^ I also have 
an asymptotically negligible lensing effect on mode /. 
If the convergence power A^Cf ~ A^C')(/(27r) is slowly 
varying on small scales A > Z, then a particular temper¬ 
ature mode I is only affected by lensing modes around 
a narrow range around A ~ y/2l. As can be seen in 
Figure n the convergence power is indeed slowly vary¬ 
ing on angular scales I > 1000. In the approximation 
that the lensing modification of the temperature power 
spectrum is small, the fractional lensing effect is roughly 
ACf/cf « Ag'‘/2, where A^'^ = s^C'”/(27r) is the 
convergence variance per logarithmic interval at scale 
s = 

The limit of the long wavelength lensing mode can be 
understood physically. Consider a slowly varying con¬ 
vergence lensing a region that is approximately scale- 
invariant on smaller scales. The lensing will result in 
slowly varying changes in scale over the region. How¬ 
ever, due to the small scale invariance, the appearance 
of the region is unchanged. 
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If the leasing mode has a wavelength somewhat smaller 
than that of a particular long wavelength temperature 
mode of the scale invariant region, the deflection field 
will interact with the fluctuation on that scale and mod¬ 
ify the power. The modification on that scale will be 
proportional to the convergence power spectrum at that 
scale. 


2.3. General Approach 

The leasing equation derived above, Eq. may 

be considered a first approximation to the full solution. 
In this section we derive the general series of differential 
equations. We have seen above that the equation govern¬ 
ing the change in the temperature correlation function is: 


ACl(Ad) = Cl+,A(A0)-Cl(A0) 

=-i(T,(/3A)T,-(/3A')>(Aci*Ad^) 


(30) 


where ,3 a = /3a(^) is the source position corresponding 
to the observed position 0, including leasing modes up 
to wavenumber A, ,3 a^ is defined analagously. For nota- 
tional convenience and clarity, derivatives of the source 
temperature with respect to the source position (3\ are 
denoted by subscripts, Ti{(3A), and derivatives of the ob¬ 
served, lensed temperature field with respect to the ob¬ 
served coordinate, 6, are denoted as diT. 

In the simplest approximation, we evaluated the right 
hand side of the equation by approximating the slopes of 
the temperature field on the unobserved source plane by 
the observed slope on the lensed image plane. 


T^if3A{e)) ^ ^T(/3a(0)) 

= Af“(0,A)AfA(0)«ftf(e) 


(31) 


where we have projected to the observed image plane and 
introduced the inverse magnification matrix M^{0,A) = 
d0 °‘which is then approximated as identity, gener¬ 
ating a relative error of leading order 0 (k). 

In order to make a better approximation, we must de¬ 
termine an equation to describe the correlator of first 
derivatives on the right hand side of Eq. (EOl)- 


{T,{[3a+sk)T,{[3'a+sa)) = mf3A + d)T,{f3'A + d')) 

= {U^a)T,{13'a + Ad)) 


= {UHa) 


TAI^'a) + 


= {TAMTA(3'a)) - -(T,,(^A)T,A(/3)v))(Ad“Ad^) 


(32) 


The change in the first derivative correlator is determined 
by the second derivatives on the source plane. A similar 
derivation can be done for the second derivative correla¬ 
tor. We get a hierarchy of coupled equations: 

AC^{A0) = -i(T,T')(Ad*Ad^') (33) 

A{T,T’) = -^{T.aT'AiAd'^AdA (34) 

A{T,aTA) = -i(T.a,rj,,)(Ad=Ad^) (35) 


and so on. The system of equations is generally not 
closed, as the change in the correlator of nth temperature 
derivatives is determined by the correlator of (n -I- l)th 
derivatives. However, at any level of the hierarchy, one 
may terminate and close the system by approximating 
derivatives on the source plane with derivatives on the 
image plane in an analogous manner as in the first ap¬ 
proximation, making an error 0{k) smaller than the 
terms retained at that level. For example, we may close 
the above system at the third equation with the approx¬ 
imation, (TiacTjhe) ~ dcd'ATiaTA)- 


2.4. Spin Decomposition 

In this section we illustrate how one can take the posi¬ 
tion space equations from the previous section and turn 
them into coupled differential equations for power spec¬ 
tra in Fourier space. The correlators of temperature 
derivatives can be broken down simply in Fourier space 
into power spectra of E and B type spin quantities. This 
approach is analogous to the spin decomposition of po¬ 
larization. First, we illustrate how to decompose the first 
derivative correlators and derive Fourier space equations. 

Define the vector v as the gradient of the source tem¬ 
perature on the source plane as a function of the ob¬ 
served position, 0, implicitly through its dependence on 
the source position /3a{6), which is also dependent upon 
the current cutoff scale A. 


vA0) = TAf3Am 


(36) 


From this the complex spin 1 and spin -1 quantities can 
be defined 




-ipl^iie 


r (fi 

v{e) = {Voo + iVy){e) = J 

i{0y = - ivy){e) = J -^^a-i{l)e 

Reality requires a-i{l) = —a*(—Z). If we define 


Eil) =-iaAl) + a.yi)) 
H(Z) = i*(ai(Z)-a_i(Z)), 
then these quantities^ can be rewritten: 


(37) 


(38) 


;( 0 ) = {vx + ivy){e) = 


dH 


[Ell)-iB{l)]e^^^e 


itj>i iie 


(39) 


(40) 

where E and B denote the F-type and H-type parts 
of the vector quantity with opposite parity. This can 
be seen by defining E{6) = J (Z^l/(27r)^ F(Z)e®* ® and 
B{G) = f d^Z/(27r)^ i3(Z)e** ®. These can be written as: 

2E(d)= J J d?0e-^'’-^ -f (6>)] 

(41) 

^ These should not be confused with the E and B modes of spin 
2 polarization quantities. 
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-2iB{e) = 


(fl 


Jio 


(42) 


Under a reflection across the x-axis, (j)i —4>h 'Vy 
—Vy, V V*, V* —>■ V. Thus we see that E is even under 
parity, E ^ E, and B is odd under parity, B —> —B. A 
similar argument may be made for reflections across the 
j/-axis. Translational and parity invariance dictate: 

{E{l)E*{l')) = {2TTfCf^5'^{l - V) (43) 


{B{l)B*{l')) = {2Tr)^CP^S^{l - I') 


{E*{l)B{l)) = 0 

From this one can form the quantities: 

{v{e)d{e')*) = (v^v^ + vyVy) 

r rfl 

- I 


(27r) 


2 I'-'l 


BBl 


(44) 

(45) 

(46) 


{v{e)v{e )) = {Va;V'^ - VyV'y + i{vyv'^ + v^v'y)) 


J {2nY 

In order to form the differential equations for the new E 
and B power spectra, one simply forms the incremental 
change (from A to A+dA) in the first derivative correlator 
{A{TiT') ) as derived in the previous section. 

A(7;r'+T,T') = [ [Cr + CP^] (48) 


. (2^)^ 

A(r,r' - Tyx; + i{TyT^ + t.t')) 
r <fi 


(27ry 


A [Cp^ - CP^] 


(49) 


Now we take the appropriate Fourier transform and take 
the limit dA ^ 0 to get equations for dCp^{A)/d\nA 
and dCP^{A)/d\nA. The initial conditions are that 
Cp^ = Pep and CP^ = 0 at A = 0, which can be 
deduced from Eq. (1^ and Eq. (|T7| by considering the 
case vpe) = T,(/3a=o(6>)) = T^e) = d,T{9). 

We now illustrate this derivation for our second order 
approximation of the system of equations: 


A{ffp = — {T,T'A{AddAdp 

A(T,rj) = -i(r,,T;,)(Ad“Ad^) 

« -^-^a^'^,{T,Tp{Ad<^Adp 

Combining Eqns. (I3ni and ISU one obtains: 

. f d‘^1 


(50) 


(51) 


(Stt)' 
1 


-{Cp^ — e** 


= -ndadhivxv'^ - VyVy + i{vyv'^ + ))(A(i“Ad'’ 


= -^dad'. 


d^l 




dAACp f d(j)A dH 


27r 


27r (27r)^ 


where we have defined for convenience Cj = Cp^—CP^. 
Taking the Fourier transform J d^Ad ' , we 

get an expression for the change in new spectra: 

d(j)A 
27r . 


A 


(cl 


EE r'BB\ _ 


C,7(UA)= 


- 


-A| 


(A • (Z'- A)) 2 e 2 *(’^‘'-A-'^<') 


(53) 


It is also convenient to define Cp^ = I "^Cp^, 
CP^ = l-^CP^, and Cf = Cp^ - CP^ so that 
Cp^ and CP^ have the same units as Cp. Note 
that Re = 2cos^((()p_A — ^n) — 1, and 

cos(0r_A - <^/) = (i - A) • - A|). Then 


\l-Ap Re 


,2i(4>i_A-4l) 


= PF{A/1, fl = cos((^A - M) 

(54) 

with F(x,ii) = 1 — — 2x^1 + 2x‘^iP. The imaginary 

part vanishes upon integration. Finally, taking the limit 
(5A —> 0 the expressions simplify to: 

-PA^Ct f d4>A 

' 2tt L' 


d In A 


27r 




(55) 


-Cp^,_^^[A/l-^^pF(A/y^J) 

A similar derivation can be done to And an equation de¬ 
scribing the evolution of Cp = Cp^ + CP^. The result 
is 

-PA^cf /-d^Ar 


——C+ = 
dlnA ' 


27r 


‘^cp 


-[A/Z-^]^G(A/Z,/r)(7+_^l 


(56) 


with G(x,ij) = \ + x^ — 2xn, and G{A/l,fj. = cos((/)a — 

^P) = \l-Ap/P. ^ 

Finally, the position-space equation for the correlation 
function, Eq. (CT . is converted to an equation for Cp (A) 
by breaking down the first derivative correlator on the 
right hand side into E and B parts. The Fourier trans¬ 
form is taken analytically and the limit dA —> 0 is taken. 
The final set of coupled equations for the second order 
solution is: 


dCP(A) 
dlnA 

-q7_A|(A) 


PA^Cp 


dcp^ 


(A) 


dlnA 


27r 

A 


PpPep 

27T 


depA 

27T - 


A“(A)d' 


.A 


Gr(A)-G(y,M)G|r^|(A) 


(57) 


dtpA 

27T . 




(G|fy^i(A)(i - 


A 

1 


-cPPaP^) 


A2 


A|CU-^ 


(58) 


dcp^ 


(A) 


dlnA 


PA'^Cp 

27r 


(1 _ giA.Ae)(^. A)2G“e2*‘^'e*''^® - 

(52) 


d(/>A 
27r . 

A 


Gr(A)M^ 


(qfyA|(A)(l-yd)^ + qt 


EE 




(59) 
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where /r = cos(0a) and v = sin((/)A). The initial condi¬ 
tions are simply Cj (A = 0) = Cj, = 0) = Cf, 

and CP^{A = 0) = 0. 

As compared to the first order equation, Eq. (P|l . 
the role of the lensed temperature power spectrum, 
Cf (A), in generating the lensing corrections on the right 
hand side has been replaced by the quantity Cf^{A) — 
{A). This spectrum generates additional corrections 
to the lensing and is determined through the second and 
third equations. Furthermore, additional corrections are 
produced through the spectrum also evolved from 

the initial condition through the auxiliary equations. It 
is through the additional dynamics of the auxiliary spec¬ 
tra as a function of the cutoff A that improvements to 
the first order solution, Eq. (P|l . are made. 

A similar decomposition can be done for second deriva¬ 
tive quantities, in order to generate the equations for 
a third-order approximation (e.g. Eqs. 123) - (|S3)- 
In this case one can decompose the second derivatives 
of the temperature on the source plane, Ty, into three 
quantities determined by the spin ±2 quantities t{6) = 

\{Txx-Tyy) + iTxy, t* { 6 ) = \{Txx - Tyy) - ITxy, und the 

scalar s{9) = ^{Txx + Tyy). In a similar fashion as be¬ 
fore, these will be described by four auxiliary power spec¬ 
tra, subject to parity invariance, C'f'^^^(A), C^'^^'^{A), 
^£ 2 ^ 2 (A), and C'f^‘^^(A) where the subscript 2 indicates 
the decomposition of second derivatives to avoid confu¬ 
sion with E and B from the first derivatives (which in 
turn should not be confused with polarization E and B), 
and S 2 denotes the Fourier amplitudes of the scalar s. 
These derivations can be done mechanically if necessary, 
in a manner similar to that given above; we will not give 
the solutions here, but mention this only to illustrate the 
hierarchical structure of this approach. 

2.5. Interpretation 

It is useful to discuss how this new formulation of 
the weak lensing calculation differs from p revi o us ap ¬ 
proache s . Th e standard approaches (e.g. IhI (Hooa), 
iCooravl (Hool) estimate the lensed temperature at a 
given observed point by expanding the unlensed tem¬ 
perature in a Taylor polynomial in the deflection angle 
(e.g. Eq. I), keeping only a few terms. The accu¬ 
racy of this expansion will depend on the smoothness of 
the underlying source temperature field, which depends 
upon how the fluctuation power is distributed over scales 
in the intrinsic temperature power spectrum. The key 
point is that the unlensed intensity is determined from 
the lensed image position in a one-step extrapolation de¬ 
pending upon the temperature derivatives evaluated at 
a single point, i.e. 6 = /3a=o(®)- 

The main difference in this new method is that the 
source position is determined in multiple steps, as a func¬ 
tion of the cutoff A. When one integrates over a new shell 
of the lensing modes with wavenumbers between A and 
A-|-(5A, the source position (3\{6) changes incrementally 
through the lens equation. The first position-space equa¬ 
tion for the temperature correlation function tells us to 
evaluate the source temperature derivatives at the new 
source position I3\ in order to determine the incremental 
modihcations to the lensed correlation function. Hence, 
by integrating over the lensing spectrum mode by mode, 


one is continuously updating the current source positions 
and the current local source gradients, which determine 
the lensing modihcations. 

This is analogous to numerical methods for extrapo¬ 
lating a function f{x) to get to f{xi) from /(xq). The 
simplest thing to do is evaluate the derivatives at the 
initial point Xq and extrapolate to get f{xi) with a Tay¬ 
lor polynomial in one-step. The accuracy will depend on 
the order of the polynomial one uses and how far away 
Xq is from Xi relative to the scale over which f(x) is 
varying. More sophisticated methods would use multiple 
smaller steps in which the derivatives of f{x) are eval¬ 
uated along the trajectory between the points xq and 
xi- Previous methods are analogous to the one-step Tay¬ 
lor expansion calculation. Our new differential approach 
is analogous to the latter, adaptive approach, in that 
the lensing modihcations are calculated along the trajec¬ 
tory on the source plane between the observed position 
/3a=o(^) = d and the hnal source temperature position 
/3a=oo(0). 

The hrst equation, Eq. ( 1201 ), which determines the 
lensed temperature correlation function as a function of 
the cutoff, is exact in the limit of the flat-sky and in the 
limit that we take the lensing shell increment (5A ^ 0. 
This is a consequence of the fact that we did not expand 
in the full deflection angle, but only in the differential in¬ 
crement of the dehection held containing the new modes. 
However, in order to evaluate the right hand side of the 
differential equation, one must evaluate the hrst deriva¬ 
tive correlator at the new source positions. The implica¬ 
tion is that we have to evaluate the second derivatives at 
the changing source position /3(A). Then to update the 
changing second derivative correlator we have to calcu¬ 
late third derivative correlators, and so on. This leads to 
the hierarchy of equations for the temperature correlator 
and temperature derivative correlators. We can termi¬ 
nate and close the hierarchy at any level by approximat¬ 
ing a derivative on the source plane with a derivative on 
the observed image plane, at the cost of neglecting terms 
representing the relative distortion of these derivatives 
through the dimensionless convergence and shear, which 
are presumed to be small (|k|, | 7 i|, I 72 I < I) in the limit 
of weak lensing. 

This brings us to a second major distinction of this 
method. Whereas the previous methods are essentially 
perturbative in the full deflection angle, so that the trun¬ 
cation errors are simply those attributed to high-order 
terms in the deflection angle, in this new formulation 
the approximation is based upon the smallness of ele¬ 
ments of the distortion tensor, which are derivatives of 
the deflection angle, in weak lensing. Intuitively, the key 
point we emphasize here is that the RMS magnitude of 
the deflection at a single point is not the physically rel¬ 
evant quantity in determining the lensing modihcation 
to the two-point function or the power spectrum. What 
is physically important is the difference between the de¬ 
flection angles at two points, which is determined by the 
dimensionless second derivatives of the projected poten¬ 
tial through the convergence and shear. Hence, in the 
context of the weak lensing calculation, it is more appro¬ 
priate to approximate about the dimensionless elements 
of the distortion tensor rather than the deflection an¬ 
gle. Furthermore, this approximation depends primarily 
on the weakness of the distortion tensor elements rather 
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than the properties of the temperature fluctuations being 
lensed. 


3. NUMERICAL IMPLEMENTATION AND RESULTS: 21 CM 
FLUCTUATIONS 

The numerical implementation of these systems of 
equations is simple. In order to implement the second 
order approximation, for example, to compute the lens¬ 
ing effect on a power spectrum sampled at N points, one 
simply creates a vector of 3iV elements. This is initialized 
to the values of the unlensed temperature spectrum, Cf 
for the first N values, the initial values of (which 
are Cf) for the next iV, and then the initial values of 
, which are zero for the final N slots. This vector is 
then fed into an adaptive Runge-Kutta 4th and 5th or¬ 
der algorithm, and an adaptive Gauss-Kronrod scheme 
is used to evaluate the one-dimensional integrals on the 
right hand sides of the equations with the help of smooth 
interpolations between the sample points. Since the in¬ 
trinsic power spectra are very smooth, these integrations 
are simple. 

We apply our first and second order sets of lensing 
differential equations to the case of 21 cm fluctuations. 
The fluctuations in the 21 cm emission at redshift z as 
seen against the CMB are described by: 


6T{v) =23 mK 


0.15 


Ts — Tcmb \ 

"Ts ) 


1 + z 
10 


0.02 


1/2 


xnil + 5) 


(60) 


(e.g. iZaldarriaga et ahl ll2004M . In the simplification 
that the spin temperature Ts is coupled to the kinetic 
temperature Tk of the gas that is heated far above the 
CMB temperature, Tg ^ Tcmb , the fluctuations are in¬ 
dependent of the spin temperature. Then the intensity 
fluctuations are due to the fluctuations in the quantity 
Ip = a:?i(l-|-(5), where Xh is the neutral hydrogen fraction, 
and S the matter density contrast. We have ignored ad¬ 
ditional sources of fluctuations, e.g. peculiar velocities, 
for simplicity. 

We have taken semi-analytic models of the expected 
21 cm fluctuation power spectra at several redshifts. 
These correspond to a sampl e model of a sing l e reion - 
ization history as described in iFurlanetto et alJ ll2004all . 
This model is described by an ionization efhciency factor 
C = 40, and we use the resulting temperature power spec¬ 
tra at redshifts of 2 : = 12, 13, 15, and 18, corresponding 
to mean neutral fractions xh = 0.25, 0.5, 0.8, and 0.96. 
These power spectra are shown in Fig. El For a lower 
ionization efficiency, reionization will occur at lower red¬ 
shifts; however, the po wer spectra remain sim ilar for a 
given neutral fraction llFurlanetto et al.ll2004ar . Hence, 
one expects that the lensing of 21 cm spectra for different 
models at different redshifts will be qualitatively similar, 
modulated by the changes in the lensing power spectra 
to those redshifts (Fig. P). 

We apply our first and second order lensing systems 
of equations to calculate the effect on the angular power 
spectrum. Solutions are shown in Figure |21 The frac¬ 
tional change in the temperature power spectrum is 
shown in Figure 01 Comparing the solutions to the first 
and second order approximations, we have found that the 



Fig. 2. — Model power spectra of 21 cm brightness temperature 
fluctuations used in this paper. Angular power spectra versus an¬ 
gular multipole are shown in the top panel. In the bottom panel, 
the three-dimensional power spectrum is plotted versus comoving 
wavenumber. The curves correspond to 2 = 12 (solid), 2 = 13 
(long dash), 2 = 15 (short dash), and 2 = 18 (dash dot). 



Fi g. 3. — Solutions to the system of differential equations Eqns. 
IS71 - EHJ as a function of the lensing cutoff A at representative 
scales. The curves describe AC/(A)/C/ = (C^^(A) — C^^)/C^^ scaled 
to the flnal value at A = A^nal = 3 x 10^, at multipole values 
I = 10^ (solid), I = 10^ (long dash), and I = 10^ (short dash) for 
the z = 12 (thin) and z = 15 (thick) model power spectra.The small 
scale modes at I ~ lO"^ where the power spectra are relatively near 
scale-invariant are comparatively unaffected by lensing modes A < 
and mostly lensed by deflection modes A > At intermediate 
and larger scales, temperature modes are affected by a window of 
lensing modes wavenumbers A ~ 10^ to lO'^. 


difference is typically of order 1% of the lensing effect. 
This demonstrates that our calculations are convergent 
and are a complete description of the lensing effect from 
degree scales to arcsecond scales. There is no evidence 
that the calculations break down at any scale, despite 
the dominance of small scale power. 

The results show that the lensing modifications are 
small on scales I > 1000 for late-time spectra {z = 
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/ 

Fig. 4. — Gravitational lensing effects on 21 cm angular power 
spectra. The fractional modifications of the angular power spectra, 
ACi/C, = {Cj-Cj)ICj are plotted as a function of angular scale 
L The curves correspond to the z = 12 (solid), z = 13 (dashed), z 
= 15 (dashed-dot) and z = 18 (dotted) models as described in the 
text. The lensing modifications towards the end of reionization 
are much smaller on small angular scales due to the flattening 
of the power spectra by the distribution of bubbles. At earlier 
epochs, the small-scale lensing is larger due to the steeper slope. 
The gravitational lensing effect on large scales is larger due to the 
steeper slopes of the power spectra. 


12,13). The intrinsic power spectra approach a scale- 
invariant slope at these small angular scales. We have 
seen that the effect of lensing on a near scale-invariant 
spectrum is small. For lensing modes that are slowly 
varying compared to a temperature mode I there is no 
effect. Since a mode in a scale-invariant spectrum is not 
affected by longer-wavelength lensing modes, the tem¬ 
perature modes at high-/ are unaffected by most of the 
lensing power at A < /; the lensing effect is determined 
by the convergence power near the scale 1. This can be 
seen in Figure 13 

Note that the largest relative lensing effect occurs 
on the rising slope of the temperature power spectrum 
PCf/{2 tt), with the most pronounced effect occurring 
on the steepest slope. The implication is that, as reion¬ 
ization occurs, and the bubble feature “shoulder” is im¬ 
printed on the angular power spectrum on scales I ^ 10^ 
to 10^, the relative lensing effect on the scale of the fea¬ 
ture and on smaller scales diminishes as the spectrum 
flattens out to near scale-invariant form. 

4. LENSING EFFECTS ON THE THREE-DIMENSIONAL 
POWER SPECTRUM 

Just as weak lensing will modify the angular power 
spectra of 21 cm fluctuations, it will also affect the ob¬ 
served three-dimensional power spectrum of the inten¬ 
sity, Pi{k). The isotropic nature of the random field 
of fluctuations at a particular epoch implies that the 
intrinsic power spectrum will be approximately spheri¬ 
cally symmetric, Pi{k) = Pi(k)', however asymmetries 
are introduced by velocity field distortions (Barkana & 
Loeb 2005). The three-dimensional Fourier tranform of 
the power spectrum is the three-dimensional autocorre¬ 
lation function of the intensities, (/(r)/(r')). ft is clear 


that the observed positions of the intensities will be dis¬ 
placed from the true source positions through the ef¬ 
fects of weak gravitational lensing. In this section, we 
will show that gravitational lensing will modify the ob¬ 
served power spectrum, and also will introduce angular 
anisotropies into P/(fe). 

We consider the scenario in which an image cube of 
intensities in a three-dimensional region is constructed 
with observations all at approximately the same epoch. 
In principle, another anisotropy will be introduced due 
to the evolution of the intensity power spectrum with 
time. Thus, the image cubes must be constructed over 
a redshift range that is narrow compared to the rapidity 
with which reionization occurs at a given epoch. 

In general the lensing power to each redshift (fre¬ 
quency) plane of the image cube will be different. The 
variation of the lensing effect between two redshift planes 
can be characterized by the convergence cross-power 
spectrum between the two redshifts, zi < Z 2 ' 

.... ^ ^ 

4c^ Jo a2(r) \x{r) J 

(61) 

The cross-correlation coefficient is then defined as: 

( 62 ) 

where ki = k(zi) and K 2 = k(z 2 ). We have calculated 
the cross-correlation for typical redshifts in this regime 
over the expected radial widths of the image cubes, on 
the order of ^ 100 Mpc. In Figure|3 the cross-correlation 
coefficient of the convergence between z = 10 and neigh¬ 
boring redshift planes up to z = 11 is shown for several 
multipoles between I — 10^ and 10®. The coefficient is 
negligibly different from unity over comoving radial dif¬ 
ferences Ar ^ 200 Mpc. The decorrelation of the lens 
convergences is only < 0.1% over this range. This is pri¬ 
marily because the lensing efficiency is greatest for lens 
planes at substantially lower redshifts of order z ^ 1, so 
that small differences in the source distance make negli¬ 
gible differences to the total lensing effect. 

Consider an image cube observed centered at red¬ 
shift z*, at a central radial distance r» and comoving 
angular diameter distance Zl*. Assume that the ra¬ 
dial width of the box is small enough so that for each 
pair of constituent planes at different frequencies (red¬ 
shifts), the cross-correlation of the lensing convergences 
is near unity. We have shown that this is an excel¬ 
lent approximation for radial widths up through hun¬ 
dreds of Mpc. The deflection field traversed by light 
from a diffuse source generally depends on its distance 
and redshift in the three-dimensional volume: S6{6;z). 
This approximation allows one to make the simplifica¬ 
tion that the deflection held seen by each of the mul¬ 
tiple frequency-planes comprising the box is the same: 
d0{9,zi) ~ S9{9 ,z2 ) « S9{9,z^) = d0{9). The statis¬ 
tics of the dehection helds can thus be described by a 
single power spectrum Cf'^{z<t). 

The apparent comoving transverse displacement due 
to lensing of an observed ray /(r) = I{9D{z),rz) 
from the plane at redshift zi is Ar±{9D{zi), zi) = 
D{zi)S9{9,rzi))- We will make the further simplihca- 
tion that the proportionality factor, the comoving angu¬ 
lar diameter distance, is approximately the same for each 


















Lensing of 21 cm Power Spectra 


13 



where we have expanded to second order in the rela¬ 
tive deflection angle. Here V_l denotes the gradient 
along the transverse direction (on the sky), and as before 

ASe = 66' - 66. With VlVi/(r) = / (Pk/{2'Kf - 

A:5_/(fe) exp(ifc • r), and defining the lensed, ob¬ 
served power spectrum, P/(fc) with {I[r)I[r'))= 
J k/ (2Tr}^ Pj(k) exp(ik ■ Ar), we find upon taking the 
three-dimensional transform, 

Pj(k^, fell) - Pj(k) = -J .0^cf{pi{k^, • If 

- Piik^ - l/D,, fc||) \l ■ (p,k^ - Z)]"} 

( 66 ) 


Ar (Mpc) 


Fi g. 5 .— Cross-correlation coefficients of the lens convergence, 
Eq. 1621 The curves represent the coefficient Ri{zi,Z 2 ) between 
2:1 = 10 and redshifts 22 between Z 2 = zi and Z 2 = 11. This range 
corresponds to a radial width of Ar = r( 2 : 2 ) — t[zi) 210 Mpc. 
The different curves represent a range of multipoles from I = 10^ 
to I = 10^. The lens convergences decorrelate only by < 0.1% over 
radial widths of 200 Mpc. 


plane in the box, and equal to the central distance. Thus, 
Ar±{6D{zi), zi) « D{zP66{6,r{zi)) ^ D^:66{6, zp. 
This approximation will be good so long as the box is 
not too long in the radial dimension. For example, in a 
box of length ~ 100 Mpc at a typical distance ~ 10^ 
Mpc, the variation of D{z) will only be of order a few 
percent across the box in the radial direction. 

The intrinsic intensity held due to huctuations within 
the box is described by 


I{r) = 


d^k 


I{k)e 


ihr 


j (27r)3 

= I{6D{z),rP = 


d^k± (Zfc|| 

(27r)2 27r 


/(fc)e""'-" e 


ikj_-dD{z) Tz 

(63) 


where we have split up the comoving coordinate r into 
components perpendicular and parallel to the line of 
sight. Here, D(z) denotes the co-moving angular diam¬ 
eter distance to redshift z, is the comoving distance 
to redshift z, and k denotes comoving wavevectors, also 
split up into transverse and parallel components. The 
ensemble average over intensity huctuations of the cor¬ 
relation function is 

/ d^k 

^Pjik)P'^ir-r) 

where Pi{k) denotes the intrinsic three-dimensional 
power spectrum. Under the simplihcations stated above, 
gravitational lensing will remap the observed intensity 
/(r) = I{6D, rP = I{D6 -|- D^66{6), rp. The observed, 
lensed correlation function is modihed as: 

(/>)/>'))7.0 = {im6 + 66),rPl{D'{6' + 66'), rp)i,^ 
= {I{D6, rpi{D'6' -P D'66' - D66, rpi^^, 

« {I{r)I{r'))i + (/(r)fviVi7(r'))/^*^(Adr 

(65) 


where in the case of spherical symmetry of the intrin¬ 
sic power spectrum, P{k_\_,kp = T’([fc^ -I- fc^]^/^). This 
equation is structurally similar to that of the hrst or¬ 
der lensing ehect on an angular power spectrum, except 
here we have a three-dimensional intrinsic power spec¬ 
trum. It is easy to see from this equation how anisotropy 
is introduced into the power spectrum through lensing; 
the right-hand side of the equation depends separately 
on the transverse component, k^. It is also apparent 
that under the approximations, the lens effect on a par¬ 
ticular mode with line-of-sight wavenumber k\\ only in¬ 
volves mixing with other modes with the same line-of- 
sight component. Thus, we may consider the case of 
three-dimensional lensing as analogous to the case of two- 
dimensional lensing on effective power spectra P{k±, kp 
as a function of k± for each fixed fc||. 

The structural similarity of this lensing formula to that 
of the two-dimensional case (Eq. (EHll) su ggests that the 
differential techniques developed in ESI can be readily 
adapted to this case in order to calculate higher order 
corrections and to verify convergence. This derivation is 
done in a nearly identical fashion. We can simply make 
the replacement (7/(A) ^ PA{k±,kp to denote the power 
spectrum lensed by lensing modes with wavenumber less 
than A, and substitute I k±D. The resulting differen¬ 
tial equation at first approximation are: 


- pA{\k^ - A/Dp,kp [A • {Dk^ - A)]'} 


(67) 


with PA=o{k±,kp = P{k±,kp. The analogous system 
at the second approximation is: 


dpAjk) 

din A 


7? l^{(D.k^^AfP^{kuk,) 


+ A^klDlP^^{kA,kp 


-ApA-D,kp^P^^i\ki_ 
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dp^^jk) _ - 
din A 

- [A • - A)]"(f-r (|fc - ;^l,fc||)(l - lArd? 


k±D^, 


+pri\k±-^m) 


A2 

'klof 


’)} 


(69) 


dPl^{k) _ PPci f d(/)A 


din A 27r 

- [A.(d4*fci-A)]"(pf^(|fc 


/ ^{(C*fcA-A)2pf^^(fcA,fc||) 




A A ^ 


)} 


(70) 


where fj, = cos((()a — (f>k_L) and v = sin(())A — (l>k±), with 

pA=o{k±,ki\) = pA^Po{k±,ki\) = P{k±,kii), P^£q = 0 as 
initial conditions. In the limit that the first approxima¬ 
tion, Eq. lHI), is sufficient, and that lensing modifica¬ 
tions are small, then the solution, Eq. ijnsji, is a suffi¬ 
cient description of the lensing effect, as one may replace 
PA(fc) ~ P{k) on the right hand side. 

In FigureElwe show the solutions to these equations for 
the models of intrinsic 21 cm power spectra as discussed 
m Here we have assumed for simplicity that the intrin¬ 
sic power spectrum is spherically symmetric and ignored 
any components, such as velocity field distortions, that 
may vary with angle cos a = /cy/fc. Plotted is the lensed 
fraction of the power spectrum, which is defined as 


A^(fe) 


Pi{k)-Pi{k) 

Pi{k) 


(71) 


as a function of the wavenumber k for several orienta¬ 
tions (dotted lines) of the wavevector k = {k^^,k±) = 
{k cos a, k sin a), where a is the angle between the line of 
sight and k. These orientations range from fully trans¬ 
verse on the sky (a = 7>'/2) to aligned with the line- 
of-sight (a = 0). In principle, this treatment is not 
appropriate for nearly line-of-sight modes, i.e. those 
with very small or no transverse component on the sky, 
k± ~ 0, because we have adopted a flat-sky formalism, 
and sky curvature effects will be important for these 
modes. In practice, these modes will not be observe- 
able due to small finite survey areas. This limit may 
be described, at each k = {kj_ -|- by the lensed 

fraction for modes with I ^ k±D = 100. This limit 
is only appreciably different from the lensed fraction of 
line-of-sight modes for modes with k ^ 0.01 h Mpc~^; 
on small scales k ^ 0.01 h Mpc“^ this limit is very close 
to the line-of-sight mode. We also compute the spher¬ 
ically averaged lensed fraction, J_^dcosa A^{k,a)/2, 
(thin dashed line), and the anisotropy ratio, which, de¬ 
fined as 


R^{k) 


1-1- max A^{k,a) 

0<oc<7t/2 


1-1- min A^{k,a) 

0<Q;<7r/2 


(72) 





Fig. 6.— Lensing effect on the three-dimensional power spec¬ 
trum P(k). The curves represent the fractional lensing effect, 
A^{k, a) = (-P(fc) — P{k))/P{k), on the power spectrum of modes 
with a given magnitude k and orientation a, as described in the 
text. These pertain to the models at z = 12 (top left), z = 13 (top 
right), z = 15 (bottom left), z = 18 (bottom right) as described in 
the text. The upper solid line is for purely line-of-sight wavevectors 
(a = 0), the lower solid line represents purely transverse wavevec¬ 
tors {a = 7r/2), and the thin dotted lines represent intermediate 
orientations at a given k, a £ (0, 7r/2). The short dashed line is the 
spherical average of the fractional lensing (over all orientations a). 
The long dashed line is a measure of the anisotrcmy due to lensing 
equal to the maximum minus the minimum of A^(fc, a) for a given 
k. 


quantifies the amplitude of the variation in lensing over 
all orientations, and measures the width of the filled re¬ 
gions in Fig. El As in the two-dimensional case, we find 
that the difference between the second and first order so¬ 
lutions is very small, and hence higher order terms are 
negligible. 

The most interesting aspect of these results is the vari¬ 
ation of the gravitational lensing effect as one rotates the 
wavevector k about a fixed axis in the plane of the sky. 
On large scales, k < 0.1 h Mpc“^, the lensing fraction 
is roughly isotropic and of order fractions of a percent. 
On smaller scales, k > 0.1 h Mpc“^, the lensed fraction 
varies significantly with polar angle a. The high incli¬ 
nation modes with small a have the largest lensed frac¬ 
tion, whereas transverse modes have near zero or slightly 
negative lensed fraction. It is not necessarily true on 
small scales that the smallest lensing effect occurs for 
purely transverse modes; this will generally depend on 
the small scale shape of the power spectrum. Meanwhile, 
the spherically averaged lensed fraction is consistently 
small, of order several tenths of a percent, over several 
decades in scale in k. The anisotropy ratio on small scales 
can approach ~l%atfc>l/i Mpc“^. 

One might expect intuitively that, since gravitational 
lensing introduces deflections transverse to the line of 
sight, only the power in transverse modes will be affected 
significantly. In fact, as is evident in Figure El the op¬ 
posite is true. The described behavior of the lensing 
effect can be understood through the shape of the power 
spectrum. In the familiar case of the lensing of the two- 
dimensional power spectrum, the quantity of interest is 
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the power spectrum of the temperature gradient on the 
sky, PCf. Analogously, in the three-dimensional case, 
what is important is the power spectrum of the tempera¬ 
ture gradient in the plane of the sky, that is klP{k^,k\\). 
In the three-dimensional case, one can consider the total 
lensing effect as the set of lensing modifications of ef¬ 
fective power spectra in the transverse wavenumber k_\_, 
each labelled by fc||. For a fixed fc||, these effective (in¬ 
trinsic) power spectra behave as 


fclP(fe) 


ik\P{k±) iffc_L3>A:||, 

\fc^P(fc||) ^ iffc_L<cA:||. 


(73) 


We see that in the case of a nearly transverse mode 
(fc_L /c||), the effective power spectrum is locally sim¬ 
ply that of the underlying power spectrum k‘^P(k). This 
is to be contrasted with the case of the highly inclined 
mode (/cj_ ^ fc||) where locally the power spectrum is 
characteristic of Poisson fluctuations: P{k±,k^^) is ap¬ 
proximately constant versus k±. The effect of this is 
to steepen the slope on scales k± fc||, thereby en¬ 
hancing the lensed fraction. We have seen in the case 
of the two-dimensional power spectrum that the Pois¬ 
son approximation ^ constant leads to lensing by 
a positive constant determined by the total convergence 
power. This result applies here approximately for highly 
inclined modes. Specifically, for k\P{k) ^ we have, 
from Eq. 

A^(fcx,fc||)«4 j = 4(«2) (74) 


so that the lensed fraction is a measure of the total power 
in convergence. For realistic spectra, this approximation 
is only valid at transverse scale k\ for lensing modes up 
to a cutoff A such that A -I- fc||, since the effective 

transverse spectrum will be shallower than Poisson at 
small scales, i.e. the large k± limit of Eq. [73 Thus, the 
above simplification serves approximately as an upper 
bound. 

Eor a given k, i.e. a vertical slice of Figure (3 the 
lensing of the nearly transverse modes can be considered 
as the lensing of the intrinsic spectrum P{k) at a high 
k± ~ k, while the lensing of the highly inclined modes 
can be considered approximately as the lensing of a Pois¬ 
son spectrum at low k±. Hence, the intuition that only 
transverse modes are lensed is incorrect: physically one 
must consider the effective transverse power spectrum 
presented to the gravitational lens. To make this clear. 
Figure | 7 | shows the effective power spectra in k± for sev¬ 
eral values of /c||. One can also see why the anisotropy 
due to lensing is small on large scales, and grows with k. 
For small fc, the effective power spectra are not apprecia¬ 
bly different as one changes the component /c||, further¬ 
more the lensing effect on the rising slope is dominated 
by contributions from smaller scale temperature modes, 
and thus relatively independent of rotating the particu¬ 
lar components (k^\,k±). On small scales or large fc, the 
k « fc|| effective spectrum is signihcantly modified com¬ 
pared to the k = k± spectrum, and the angular range, 
0 > a > 7r/2, spans a wider difference in k±. 


4.1. Relation between 3D and 2D Lensing 

One may observe that the lensed fraction of the three- 
dimensional power spectrum can become relatively large 


CL 

coH 



0.01 0.10 1.00 10.00 0.01 
kx (h Mpc'^) 


0.10 1.00 10.00 
k|| (h Mpc'') 


Fig. 7. — Left: Effective transverse power spectra in kj_ 
for the z = 12 model. The solid line is simply the three- 
dimensional spectrum /c^P(fc)/(27r^), or equivalently the effective 
transverse spectrum for fc|| = 0, e.g. fully transverse modes. 
The dashed lines are the effective spectra for (left to right) A;|| = 
0.1, 1.0, and 10 h Mpc~^. The normalization is arbitrary. For 
highly inclined modes, i.e. kj_ A:||, the effective spectrum that is 
subject to lensing is steepened compared to the three-dimensional 
spectrum (solid line) into a Poisson form. Right: Weighting func¬ 
tion, p{k±,k\\) as a function of k\\. Because of the decline in the 
power spectrum P{k) on small scales, the lensing of the inclined 
(q; > 7r/4 or fc|| > fcj^) three-dimensional modes are discounted 
relative to transverse modes in generating the lensing effect on the 
angular power spectrum. 


(> 1%) for some modes at small scales, while the two- 
dimensional lensing at equivalent angular scales is com¬ 
paratively small. To understand this, we may relate 
the lensing effects on the three-dimensional brightness 
distribution to that of the two-dimensional power spec¬ 
trum. In the limit of the flat sky"^, the two-dimensional 
power spectrum at a specific redshift {C'^{A6) = 

J d^Z/(27r)^ Cf is related to the three-dimensional 

spectrum {{T{r)T{r')) = j d^k/{2-Kf P{k)e^'^-^^) as 
Cf = (ifc||/(7rU^) P(^/U, fc||). That is, for a pro¬ 

jected wavenumber I, the two-dimensional spectrum is 
simply the sum of the power over all A:|| for all modes 
with transverse component k^ = l/D. From this a 
relation between the three-dimensional lensed fraction, 
A^(fc), and the two-dimensional lensed fraction, defined 
as Ap = {Cf — Cj)ICj, can be derived: 

poo 

Af= dfc||A^(Z/Afc||)p(V^,fc||) (75) 

Jo 


P{l/D,k\\) 


Pjl/D.kii) 

/P°*||P(l/Afc||) 


(76) 


where p{k±,k\\) is a weighting function proportional to 
the three-dimensional power spectrum and normalized by 
the total variance in modes with transverse component 
k± = l/D, such that dk\\p{l/D,k\\) = 1. 


The corres£onding_ex£ressi£naccounting for all-sky curvature 
is given in|ZaMarriuga^_et_aJJ 1200J). 
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From this equation, one can understand how the two- 
dimensional lensing effects are generated from the three- 
dimensional lensing distortions. We have seen how 
the three-dimensional lensing involves the mixing of the 
power in different k± modes for each group of modes with 
the same fixed A;||. The two-dimensional lensing effect is 
then simply the weighted sum of the lensed fractions of 
different A:|| modes with the same fixed k± = l/D. In Fig¬ 
ure | 7 | we show the normalized weight function p(k^,k\\) 
for several values of fc||. The shape can be understood in 
the same way as the effective lensed spectra are under¬ 
stood, but reversing the roles of fcj_ and k^\. On scales 
k\\ <C k±, the weight is p ~ constant with respect to fc||, 
while on scales fc|| ^ fcj_, p ~ P(A:||), which is a declin¬ 
ing function of wavenumber on small scales. At a given 
transverse scale k± = l/D, the weight function p clearly 
favors contributions from modes Qf > 7r/4, and discounts 
more inclined modes a < 7r/4. Thus, the magnitude of 
the two-dimensional lensed fraction of Cf is primarly an 
average of the three-dimensional lensed fractions for the 
near transverse wavevectors, and is not affected by the 
highly inclined modes that are more strongly modified 
through lensing. Finally, note that the above formalism 
will be inaccurate at scales I < 100 because of sky cur¬ 
vature effects, however the intuition will be the same. 

4.2. Angular Structure 

In this section, we examine more closely the angular 
structure of the three-dimensional lensed fraction. In 
Figure IHl we show the lensed fraction as a function 
of the orientation angle a for two values of k. There is no 
obvious analytic dependence on the angle from the solu¬ 
tion, Eq. However, there is a fairly smooth transi¬ 

tion of the lensed fraction from highly inclined modes to 
nearly transverse modes. This suggests a simple expan¬ 
sion in a few Legendre polynomials Pi{p) in p = cos a. 

A^(fc,a) = Ao(fc)+A2(fc)P2(p)+A4(fc)P4(p) (77) 

Here we have expanded only in even polynomials since 
the lensed fraction is symmetric around a = 7r/2, and 
the coefficient functions A 2 „(fc) are determined numeri¬ 
cally using the completeness relations. Shown in Figure 
Oare the separate components, A 2 n, and their sum. It 
is evident that this truncated expansion well describes 
the angular dependence of the lensed fraction. This ex¬ 
pansion is equivalent to a fourth order polynomial in p. 

A^(k, a) = Boik) + H2(fc)p2 + B4k)p^ (78) 

with Bo{k) = Ao(fc) — (l/2)A2(fc) -I- (3/8)A4(fc), i? 2 (fc) = 
(3/2)A2(fc) - (30/8)A4(fc), and B^k) = (35/8)A4(fc). 
T hese coefficients a r e plot ted for an example in Figure 0 

iBarkana fc LoebI ll200f)H have shown that accounting 
for redshift space distortions due to peculiar velocities 
leads to angular anisotropic components in the intrinsic 
power spectrum that can be expressed as a polynomial 
in p = cos a, i.e. 

P(fc) =P^o(fc)+P^2(fc)p2 + P^4(fc)p4. (79) 

The measurement of the separate components allows the 
separation of the various astrophysical effects that con¬ 
tribute to the 21 cm fluctuations. We have shown that 
gravitational lensing of the purely isotropic component 



Fig. 8. — Angular dependence of the three-dimensional lensed 
fraction. Here we illustrate the numerical solution (crosses), the 
three Legendre polynomial components, P2n(M = cos a), and 
their sum (solid line). This pertains to the z = lb model for 
fc = 1 h Mpc“^ (top) and k = 10 h Mpc“^ (bottom). 



Fig. 9.— The angular components of the lensed fraction of the 
three-dimensional power spectrum. This pertains to the z = lb 
model. The coefRcents are defined such that the lensed fraction is 
A^(fc) = Ao{k) + A2{k)P2{u) + Aiik)P4iiJ.) = Boik) + B2ik)u? + 
Biik)/!'^, where fi = cos a. 

{P{k) = P^o{ky) generates components with the same 
angular dependences, with magnitudes described by the 
coefficients Bn{k). This effect arises from the different 
effective transverse spectra fc^P(fc_L,/c||) that are lensed 
for different orientations of k. Future measurements of 
the p” components of the power spectrum will need to 
take into account the gravitational lensing contributions 
to each component when precisions better than half a 
percent are achieveable for scales fc < 1 h Mpc“^ and 
when precisions of order one percent are achieveable for 
scales A: > 1 h Mpc“^. 

4.3. Lensing of Intrinsic Anisotropies 

In the previous sections we have described the struc¬ 
ture of weak lensing effects on a spherically symmetric 
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power spectrum. This approximation will be most rele¬ 
vant late in reionization when fluctuations from reionized 
bubbles will dominate over angular anistropies caused by 
peculiar velocity distortions. Early in reionization, when 
bubbles are less important, the angular components that 
vary with and can be of the same order and larger 
than the spherically symmetric part. In this case, in 
order to describe the gravitational lensing effects, one 
must consider the full anisotropic intrinsic power spec¬ 
trum, Eq. (P|l. Although on the surface this appears 
much more complicated, the equations derived in previ¬ 
ous sections to calculate the lensing are straightfowardly 
applicable to any three-dimensional power spectrum with 
a general P{k±,k^^) dependence. 

We have considered lensing of an early 21 cm model, 
in the case where P^ 2 {k), P^i(k) > P^o{k). We Hnd 
that the lensed fraction is of the same order of magni¬ 
tude as those described in the spherically symmetric case 
(Fig. EJ. The overall effect is smaller than one percent 
on scales A: < 1 h Mpc~ , but highly inclined modes can 
be lensed > 1% on scales A: > 1 h Mpc~^. The lensing 
again generates angular components depending on 
with amplitudes that grow with increasing k. The an¬ 
gular profiles are qualitatively similar to those shown in 
Figure |H1 however we find that higher order Legendre 
polynomials, i.e. P;(/r) with I > 4, are needed to fully 
describe the angular variation of the lensing at smaller 
scales. 

These characteristics can be understood qualitatively. 
We expect that the order of magnitude of the lensing 
to remain the same as long as the P ^2 (A;) and P^i {k) 
components are also smooth and do not add sharp fea¬ 
tures to the overall spectrum. Next, we note that the 
lensing in the regime of highly inclined modes at large 
k\\ ^ A:_l is the same as described in the spherically 
symmetric case. This is simply because, in this regime, 
= cos^ a = A:|/(A:^ + ^1) ~ so that the intrinsic 
power spectrum is P{k) « P^o(A:||) -I- P^ 2 (A:||) -|- P^4(A:||), 
since k = {k\ + « k\\. Thus, the transverse gra¬ 

dient spectrum, k\P{k) ^ k'j_ is effectively Poisson for 
k± A:||, and hence will be lensed by a positive constant 
determined by the convergence as before. This effect be¬ 
comes most pronounced for very large A:||. 

The more complicated structure of the lensing is in¬ 
tuitively a reflection of the more complex intrinsic spec¬ 
trum. Just as the lensing of the purely isotropic com¬ 
ponent generates lensed components varying as and 
/r'^, we may expect that lensing an intrinsic spectrum 
with significant and pieces will generate higher 
order angular moments. The general structure of the 
lensed fraction is determined by the shape of the trans¬ 
verse gradient spectrum. For a highly anisotropic intrin¬ 
sic spectrum, there are three distinct regimes that de¬ 
termine the overall shape of A:5^P(A:j_, A:j|) with respect 
to A:_l, for a given slice of modes with hxed A;||. In 
the first regime k± <C A:||, and the spectrum is locally 
Poisson. In the second regime, k± ^ A:||, the factor /i^ 
is beginning to decrease but is still a fraction of order 
unity. Here, the spectrum is determined by the shape 
of the k\{P^ 2 {k)ii'^ -\- P^4{k)^‘^) components versus k±. 
In the third regime, k± ^ A:||, so that the terms 
are tiny and suppress the contributions from P ^2 (A;) and 


P^4{k). Here, the modes are mostly transverse on the 
sky, and the transverse gradient spectrum is dominated 
by the isotropic term, k\P{k) k\P^o{k±_). These 
three regimes and the transitions between them account 
for the additional structure that can be seen in the lensed 
fraction. 

5. DETECTABILITY 

In this section we discuss the magnitude of the grav¬ 
itational lensing effects on the three-dimensional power 
spectra in relation to the expected sensitivites of planned 
radio array experiments. We have seen how the lensing 
effects are typically < 1% on scales A: < I h Mpc~^. The 
first generation of experiments is unlikely to achieve this 
precision in the near future, and thus lensing effects will 
be unimportant. We consider here the potential sensi¬ 
tivity of a SKA experiment to an EoR signal at z = 8, 
with a mean ionized fraction of 0.54, and we only con¬ 
sider the lensing of the spherically symmetric part of the 
power spectrum P{k). The parameters of the SKA are 
not fixed; here we assume an array of 5000 antennae with 
a total effective area of 6 x 10® m^, and a minimum fre¬ 
quency resolution of Az/ = 0.01 MHz. The density of 
antennae is assumed to fall off as r~^ from the center. 
The total bandwidth of the observation is 6 MHz and 
the integration time is 1000 h ours. The system tem per- 
ature is set to T^ys = 440 K llBowman et al.ll^005|l and 
the effective survey size is 7r5.6^ deg^. 

In order to compute the potential sensitivity to the 
lensed three-dimensional power s pectrum P(k), w e fol- 
low the methods de scribed in iMorales) ll2005fl and 
iMcQuinn et al.l ll2005ll . The contributions from sample 
variance and thermal noise in visibilities to the uncer¬ 
tainty in the power spectrum dP{k, a) in annuli at fixed 
k and a are computed in logarthmically spaced bins in 
k. For simplicity we have ignored here the complicated 
matter of foregrounds. Next the uncertainty in the com¬ 
ponents of the model P{k) = Xi{k) + X 2 {k)^^ + Xz{k)ii^ 
is computed. Here Ai(A:) = (I -|- Bo{k))P{k), X 2 {k) = 
B 2 {k)P{k), and Xz{k) = Bi{k)P{k). The Fisher matrix 
is 

_^ 2^ Qp Qp 

F^3 ^ 5P{k,af 'dM^ 

An estimate of the uncertainty on the components Xi is 
In Figure Uni these estimates are plotted. The 
component contributions from gravitational lensing are 
smaller than the expected sensitivities by at least a fac¬ 
tor of 10 on scales k ~ 0.3 Mpc“^. If future experiments 
are to detect the gravitational lensing modihcations of 
the 21 cm power spectrum, they will need greater angu¬ 
lar and frequency resolution in order to probe the regime 
A: > 1 h Mpc“^, where lensing effects become relatively 
larger and lensing contributions to the angular compo¬ 
nents become more pronounced. If futuristic experiments 
are able to probe the power spectrum components to 
< 1 % precision, gravitational lensing effects will be im¬ 
portant to take into account. 

6. SUMMARY AND CONCLUSION 

In this paper, we have addressed the effects of grav¬ 
itational lensing of large-scale structure on the bright¬ 
ness temperature power spectrum of high redshift 21 
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Fig. 10. — Sensitivity estimates for a 21 cm power spectrum at 
z = S. The very thick solid line is the isotropic component of the 
power spectrum P{k). The thick solid, dashed, and dash-dotted 
curves labelled /r®, ij?, and correspond to the contributions 
due to lensing, Bci[k)P[k), B 2 {k)P(k), Bi{k)P(k), respectively. 
The thin solid, dashed, and dash-dotted curves are the sensitivity 
estimates for each component respectively. Here we have assumed 
an SKA configuration as described in the text. 


cm. Gravitational lensing presents a secondary source 
of fluctuations caused by density perturbations between 
emission and observation, and reprocesses intrinsic fluc¬ 
tuations caused by density, neutral fraction, and spin- 
temperature fluctuations in the EoR. 

We have presented a new formulation for calculat¬ 
ing the gravitational lens effect on temperature power 
spectra of diffuse backgrounds. This harmonic-space ap¬ 
proach leads to a hierarchical system of equations that 
describe the lensing effect on the power spectrum as one 
integrates differential shells of lensing modes regulated 
by a single cutoff scale A. We have described the general 
structure of the system of equations, which can be closed 
at the chosen level by approximating the magnification 
matrix. This method is an alternative to previous series 
solutions, which require evaluation of high-dimensional 
integrals or multiple transforms between position space 
and harmonic space. 

We have applied this method to the case of 21 cm fluc¬ 
tuations from the era of reionization. The effect on two- 
dimensional spectra Ci at a single frequency (redshift) 
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is typically less than 1% on scales 100 < I < 10®. In 
particular we And that higher-order corrections are neg¬ 
ligible on all scales despite the dominance of small scale 
power. At late times when the spectra are dominated by 
the ionized bubbles, the flattening of the power spectrum 
on small scales suppresses the lensing effect due to near 
scale-invariance of the spectrum. 

For small surveys, the three-dimensional power spec¬ 
trum P{k) can be directly probed by mapping frequency 
to line-of-sight distance. We have shown how the ex¬ 
pected lensing effect on the three-dimensional power 
spectrum can be reduced to a two-dimensional calcula¬ 
tion for each set of temperature modes with the same 
line-of-sight component wavenumber A:||. Since the rele¬ 
vant quantity is the transverse gradient power spectrum, 
which is different for each A:||, gravitational lensing will af¬ 
fect transverse wavevectors differently from line-of-sight 
wavevectors. The result is a dependence of the lensed 
fraction on the angle between k and the line-of-sight that 
is well described by a quartic polynomial in the cosine. 
Thus, weak lensing generates angular components of the 
power spectrum with amplitudes of order one percent or 
more on small scales, fc > 1 h Mpc~^. 

Finally, we have compared the predicted lensing sig¬ 
natures to the potential sensitivity limits of SKA obser¬ 
vations of the 21 cm emission. Future experiments will 
likely need greater sensitivity and resolution in angular 
and frequency space in order to detect the effects of weak 
lensing. The gravitational lensing effects on power spec¬ 
tra will become important for futuristic high-precision 
21 cm experiments that can measure the power spec¬ 
trum to better than one percent and can probe scales 
fc > 1 h Mpc“^. 
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